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_ Abstract. Adapting the recently developed randomized dyadic structures, we introduce the 

, notion of spline function in geometrically doubling quasi-metric spaces. Such functions have 

interpolation and reproducing properties as the linear splines in Euclidean spaces. They also 
have Holder regularity. This is used to build an orthonormal basis of Holder-continuous wavelets 
with exponential decay in any space of homogeneous type. As in the classical theory, wavelet 
, bases provide a universal Calderon reproducing formula to study and develop function space 

• theory and singular integrals. We discuss the examples of spaces, BMO and apply this to 
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a proof of the T(l) theorem. As no extra condition (like 'reverse doubling', 'small boundary' 
of balls, etc.) on the space of homogeneous type is required, our results extend a long line of 
works on the subject. 



1. Introduction 



The main goals of this paper are two-fold: the construction of orthonormal, Holder-continuous 
wavelet bases in general spaces of homogeneous type, and their applications in the theory of 
singular integrals and function spaces in the same general set-up. Despite several existing results 
CO ■ of related nature, the scope of our theory is completely new in at least two respects: 

^ . First, as far as we are aware, it seems that we offer the first construction of an orthonormal 

wavelet basis, as opposed to a frame, in this setting. Only relatively recently, Deng and Han wrote 
1^ . [TT| p. 40]: "Orthonormal wavelet bases are out of reach on a space of homogeneous type. Instead 

ly^ ■ the theory of frames will be used. Roughly speaking using a frame means that you tolerate a limited 

amount of redundancy while redundancy is completely avoided with a basis." Here, we build a 
genuine basis, and so avoid this redundancy. Our construction starts from appropriate splines, 
which also seem to be new on an abstract quasi-metric, even metric, space, and of independent 
interest. 

Second, we are careful not to impose any additional assumptions other than those defining a 
, , space of homogeneous type in the sense of Coifman and Weiss [S] — that is, a set X equipped with 

r> , a quasi-distance d and a Borel measure jj, that is doubling on the quasi-metric balls — , and we 

■ keep working with the original given d, without changing to an 'equivalent' one (see below). The 

difficulty is that the quasi-distance, in contrast to a distance, may not be Holder-regular, and 
quasi-metric balls might not be open, nor even Borel sets with respect to the topology defined by 
the quasi-distance (openness is assumed in |6]). We merely assume that balls be Borel sets, but 
even this assumption can be easily remedied by. See Section [S] 

In contrast to this, we emphasize that a substantial part of the rich literature on analysis in 
spaces of homogeneous type involves additional assumptions on the underlying space. Let us con- 
sider in particular the existence of a Littlewood-Paley decomposition (or a Calderon reproducing 
formula) with a regular kernel, which has been the object of numerous works. (Decompositions 
with discontinuous kernels arising from martingale differences are well-known even in more general 
set-ups; this is a different story.) Indeed, the regular Littlewood-Paley decomposition is the core 
formula that allows the various characterizations of Hardy spaces, the development of function 
spaces of Besov and Triebel-Lizorkin type, and the analysis of operators that act on them. (Here, 
we have in mind a definition of Besov spaces linked with the modulus of continuity and thereby 
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the (quasi-)metric of the underlying space; some other generaHzations from the point-of-view of 
abstract approximation spaces, as in |26) . are yet another story.) 

As far as we understand, there has always been some kind of further hypothesis either on X, d 
or /X. The seminal work in this context is the paper of David- Journe-Semmes on the T(b) theorem 
[9]. The starting point of their theory is the existence of a new Holder-continuous quasi-metric 
adapted to the measure fi and topologically (but not geometrically) equivalent to the original d, 
which is provided by work of Macias-Segovia |25l Theorems 2 and 3] (for more precise versions 
of the Macias-Segovia results, see [M] and [3T]). Their singular integrals are then defined relative 
to this new quasi-metric, which leads to a possibly different class of Calderon-Zygmund operators 
than those defined in terms of the original d and /j,. (Even then, [3] still needs some more technical 
assumptions such as no point masses, infinite mass for X and even a "small boundary" property for 
balls.) Working with this kind of new quasi-metric has become a common set-up in the literature. 
The theory of function spaces under this assumption has been developed by Han-Sawyer |16j . see 
also the nice review by Han- Weiss [17) . 

By [25], Theorem 2 (for an elementary proof, see [34]), it is possible to only change the quasi- 
metric to a metrically equivalent one that is Holder-regular. This does not change the classes of 
singular integral operators (up to changing the Holder-exponents), but leaves another difficulty: 
some estimates require a control from below for the growth of balls. See the discussion below 
about (|l.ip . In several more recent papers [T31 [HI HH [13] on function spaces on (quasi-)metric 
spaces, this issue is circumvented by assuming the 'reverse doubling' property 

n{B{x,Cr)) > {l + e)ijiB{x,r)), < r < diam(X)/C, 

which is equivalent to the non-emptyness of annuli: B{x, R) \ B{x, r) ^ whenever R/r is 
large enough and r is at most a fraction of diam(X). This latter condition is another common 
assumption; it was only recently eliminated from some results about positive integral operators 
on metric spaces by Kairema |22| . 

The reverse doubling excludes in particular the presence of point masses (which is also frequently 
assumed, even if empty annuli are otherwise tolerated) , and therefore rules out some basic examples 
like discrete groups Z or Z/pZ and the multidimensional analogs, the typical discrete metric 
structures arising in theoretical computer science (trees, graphs, or strings from a finite alphabet), 
Qp from arithmetics, or discrete approximations of other spaces of homogeneous type (like those 
constructed in |2]), even if the original space did satisfy the reverse doubling property and some 
ad hoc orthonormal (or bi-orthogonal) wavelet bases had been constructed. To mention a few 
references in these directions, see [3], [3S], [37] ■ We remark that even if Holder-regularity on 
discrete structures might not be an issue at first glance, it becomes one thinking the discrete 
structure as approximating a continuous one. 

In contrast, we repeat, we do not make any such additional assumptions or changes to a different 
quasi-metric. 

Let us come to the construction of spline wavelets. Splines have a long history in approximation 
theory, due in particular to the interpolation property and their polynomial nature, and were 
popularized by early books in the I960's. (The first book listed from a MathSciNet title search 
"spline" is [T] and the reviewer dates the appearance of the splines to Schoenberg's work during 
World War II; see [36].) As for spline wavelets, they were first constructed on the real line and on 
the 1-torus by Stromberg [38]. They were rediscovered independently by Battle [4] and Lemarie 
[53] in Euclidean spaces. Using the Multiresolution Analysis scheme of Mallat [57] and Meyer |28j . 
it became clear that the existence of regular, compactly supported splines with the interpolation 
property and reproducing formula leads to bases of regular spline wavelets. This was attempted 
in other geometrical contexts. For example, splines and spline wavelets were constructed in |24| 
on a stratified Lie group following a minimisation procedure using the existence of a sublaplacian, 
and in [7] on compact manifolds like spheres. We do not claim to be exhaustive in this history. 
Nevertheless, we have to follow a new route as there is no group structure attached, nor any local 
coordinates to help. We rather use a probabilistic approach as described next. 

Central to our analysis are the dyadic structures in a geometrically doubling quasi-metric space 
as constructed by Christ [5] (a construction in Alilfors-David spaces was done earlier by David 
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[S]), and their randomized versions recently studied by the second author with Martikainen |21| 
and Kairema |TS]. (Another variant is due to Nazarov, Reznikov and Volberg [32].) We build the 
splines as averages of the indicators of dyadic cubes under a random choice of the dyadic system. 
This can be motivated by an appropriate point of view at the classical piecewise linear splines on 
M, which are generated by the function 

s{x) = a:l(o,i](x) + (2 - x)l(i_2) (x). 

Observe that random dyadic intervals of sidelength 1, in the sense of Nazarov, Treil and Volberg 
|33|, Sec. 9.1], can be defined by translating the standard intervals [fc, fc + 1) by a random number 
u £ [0, 1). Thus the unit cube with left end at the origin, [0, 1), is translated to [u,u + 1). It is 
easy to check that the average of its indicator function i[u.u+i)i^) = l[o.i)(-'^ ^ ^) over the uniform 
random choice of u is precisely 

-1 

1[04)(^ - u) du = l[o,i) * l[o,i)(-'i^) = s(x). 

So one can indeed think of splines as averages of the indicators of random dyadic intervals. This 
is the basic idea which will guide us in constructing splines for a quasi- metric space X. 

Once the splines are available, and only at this point we add into our considerations the doubling 
Borel measure, we can obtain the orthonormal bases of scaling functions and the wavelets by the 
general procedure of Meyer. These will be Holder-continuous functions -0^ with an exponential 
localization in the ^'""-neighbourhood of a given point y^: 

for some a G (0, 1] that depends only on the quasi-triangle constant (with a = 1 in the metric 
case or when the quasi-distance is Lipschitz-continuous). The construction of boundedly supported 
wavelets in this generality remains an open problem. 

There is an important feature related to the collections of points ^V'^ := {i/q}, which index the 
wavelets of the scale S''. These collections arise as = ^"'=+1 \ where • • • C JT*^ C C 
• • • is an increasing sequence of point-sets, and each J^T*^ is both (5'^-separated and (up to constant) 
^'^-dense in X. The important feature is that '3^'^ can be much sparser than ^'^'^^ in some regions 
of the space, which reflects the 'absence of the length scale 5'"'+^' in the local geometry of the 
space. Accordingly, the distance d{x, '3^^) of a given element x € X io the nearest point S 
will be a significant quantity, which can in general be much larger that 5^ . 

This quantity appears in central technical estimates when bounding series of the following type, 
which naturally arise in the context of function spaces: 

1 _ ( __ (d{x,W^)Y\ ^ C 



^LiBix,Sf)) ' \ 'V II- t4Bix,r)) 

This (non-obvious but valid) estimate serves as a replacement of the following bound, which is 
repeatedly applied by Han, Miiller and Yang in the reverse doubling context (cf. [15[ Lemma 3.5]): 

V- 1 C 



Indeed, this is a quick consequence of reverse doubling, but invalid in general spaces of homoge- 
neous type. A typical cause to destroy this latter estimate is the presence of a large empty annulus 
around x, but this then results in some large values of d{x, 3^-'), which compensates for the failure 
of Ijl{B{x,5^)) to grow fast enough in (jl.ip . 

Although we do not develop the theory of function spaces in the detail of [T5] here, it is clear 
that our analysis involving the quantity d{x, '3^^) will extend to further related questions that we 
do not explicitly deal with; we believe that it could be used to extend large parts of the recent 
theory of 'RD (reverse doubling) spaces' to general spaces of homogeneous type. As an illustration, 
we give a short treatment the r(l) theorem using our construction. 
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Over the last decade or so, it has been discovered that several aspects of harmonic analysis 
can even be pushed beyond the setting of doubling measures; see e.g. [551 HU[ HT] for some of the 
pioneering developments in the context of E^, and [101 [5T] for recent developments in abstract 
quasi-metric spaces. We do not address the non-doubling measures here, as it seems that our 
regular splines and wavelets are most useful in the (already quite general) doubling situation. In 
fact, while our splines are constructed completely independently of any underlying measure, it 
turns out that they automatically have a good theory with respect to any doubling measure 
on the space. This would not be the case for a non-doubling measure, and a successful wavelet 
theory for such measures, if possible at all, should probably be adapted to the particular measure 
in a more complicated manner, already at the early point, where we can manage with a purely 
geometric construction. 

About notation. We use C to denote positive constants, whose value may change from one 
occurrence to the next. We also abbreviate '< C x • • • ' to just '<'. We use 7 in a similar role as a 
positive exponent. In contrast to C . it typically decreases from one occurrence to the next. For the 
measure ('volume') of balls B{x^ r) := {y € X] d{x, y) < t"}, we sometimes use the abbreviations 

V{x, r) := fiiB{x, r)), F(x, y) := V{x, d{x, y)). 

Acknowledgements. This work started during a visit of the second author to Paris-Sud whose 
support is gratefully acknowledged. The second author was also partially supported by the Acad- 
emy of Finland, grants 130166, 133264 and 218148. We thank Detlef Miiller for answering our 
questions on his work on function spaces on Reverse Doubling spaces, and the anonymous referee 
for comments to improve the presentation. 

2. Preliminaries on dyadic systems 

Introduction and motivation. In this section we review the relevant-for-us parts of M. Christ's 
[S] construction of dyadic cubes in a space of homogeneous type, as well as its recent probabilistic 
version from [TO1|21| . For the convenience of the reader, and since we feel that we have managed to 
slightly streamline the earlier presentations, we will give a self-contained treatment, even though 
the actual novelty is only in certain details of the construction. 

Christ's cubes Qjj (where k designates the generation or length scale of the cube, and a indexes 
the cubes inside a given generation) are defined by two sets of auxiliary objects: the centre points 
x^, which determine the rough position of individual cubes, and a partial order < on the family 
of the index pairs (fc, a), which determines the set inclusions among different cubes. 

In [T9l [2T] , a systematic method of constructing several Christ-type families of dyadic cubes 
Q^{uj) was introduced, where lo belongs to a parameter space fi. Here can be equipped with 
a probability measure, which gives rise to a notion of random dyadic cubes. This allows us to 
compute averages over a random choice of the cubes, and make probabilistic statements about 
them. As usual, averaging has a smoothing effect, and the random cubes will 'on average' enjoy 
better regularity properties than any fixed cubes would do. This is essential for our construction 
of Holder-continuous splines. 

Since the cubes are determined by the dyadic points x^ and the partial order <, the 
construction of the parametrized family of Q*^{uj) amounts to defining appropriate parameter- 
dependent points Zq(w) and a parameter-dependent partial order This is achieved by first 
fixing the reference objects x^ and < as in Christ's original work, and constructing z^{ijj) and <^ 
as their parametrized perturbations. 

All these constructions involve certain arbitrary choices, which all work equally well for the 
applications considered in the earlier papers. A feature of our presentation here is that we are 
going to insist on somewhat more specific, rather than arbitrary, choices in certain details of the 
construction, as this will be a convenience later when using the dyadic cubes in the building of 
our splines. 

We now turn to the details. 
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General assumptions. In what follows, X is a set equipped with a quasi-distance with quasi- 
triangle constant Aq > 1, namely, d{x,y) = d{y,x) > , d{x,y) = if and only if x = y, and d 
satisfies the quasi-triangle inequality 

d{x, y) < Ao{d{x, z) + d{z, y)). 

We assume that X has the geometric doubling property^ namely that there exists a natural number 
N such that any given ball contains no more than TV points at quasi-distance exceeding half its 
radius. This is a certain finite-dimensionality requirement on the space: for example, the infinite- 
dimensional Hilbert space £^ fails this property, since all the half-unit vectors ^Cfe, k S N, belong 
to the unit-ball, while their mutual distance is > 1/2. We do not need the doubling measure, 

or in fact any underlying measure, at this stage. 

Let us fix a small parameter (5 > 0. For example, it suffices to take 5 < yjJ^q^o Roughly 
speaking, the point is that phenomena on the length scale S^'^^ should remain much smaller than 
the length scale 6^ ^ even after repeated use of the quasi-triangle inequality where we 'lose' the 
constant Aq at every application. For example, if the points xq, a;i, X2 ■ ■ ■ satisfy d{xi_i, x.i) < (5'^'+^, 
we can still conclude that d{xQ,Xr) is much less that 5'' for r < 10. 

Reference dyadic points. For every fc e Z, we choose a set of reference dyadic points x'^ as 
follows: For fc = 0, let ,S^'^ := {x°}q. be a maximal collection of 1-separated points. Inductively, 
for k e Z+, let jr'= := {a;^}„ D ^'^-i and Jf-'^ := {x-''}^ C ^"('=-1) be maximal (5*=- and 
(J^'^-separated collections in X and in respectively. 

2.1. Lemma. For all k ^ "L and x ^ X , the reference dyadic points satisfy 

d{xl,x)i) > 5^ {a ^ 13), d{x, SC^) = mmd(2;,2;^) < 2Ao<5^ 

Proof. The separation property is part of the construction. By maximality, it follows that for all 
X <^X and fc > 0, 

JT'^) = mind(a;,a;^) < . 

OL 

Also, given x G X, we can recursively find points x^^ , x~^^, . . . , such that 

d(x,a;°J<l, d{xl^,x~\)<b-\ ... d{x-t~'-\x'^^)<b-\ 

and hence 

d{x,x'J:) < ^ < < 2A,b~K □ 

i=o ° 

Note that J^T*^ C ^'''+-'^, so that every x\ is also a point of the form a;^^^, and thus of all the 
finer levels. We denote 'S/^ := 3^^^^^ \ JT'^, and relabel these points as = {y^la. These points 
will play an important role as a parameter set of our wavelets, to be constructed. 

Reference partial order. We set up a partial order < among the pairs (fc, a) as follows: Each 
(fc + 1, f3) satisfies (fc + 1, /3) < (fc, a) for exactly one (fc, a), in such a way that 

(2.2) d(x^+i,x^)<-l-<5^- =^ (fc + l,/3)< (fc,a) =^ d(x^+\ x^;) < 2Ao5^ 

The pairs (fc + 1,13) with (fc + < {k,a) are called the children of (fc,a). Geometric doubling 
implies that their number is uniformly bounded. So far, we have essentially followed the original 
construction of M. Christ [5], a slight nuance being the choice of the point-sets in such a way 
as to have the nestedness C ^''+^, which was not required in [5l [T9l [2T]. 
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Labels for the points. For a successful perturbation argument to construct tlie parametrized 
dyadic points z^{uj) below, we need certain book-keeping among the near-by dyadic points a;^ of 
the same generation. 

Points {k,a) and (fc, /3) are called neighbours^ if they have children {k -f- 1,7) < (fc,a) and 
(k + 1,7]) < {k,l3) such that d{x^+'^ ,x^+^) < {2Aq)-'^5''. In this case 

4) < Aod{xi,x';+') + Ald{x';+\x':,+') + Aldix':+\x'^p) 

< 2Al5^ + \Aq5'' + 2Al5^ < bA^S''. 

The number of neighbours that any point can have is also uniformly bounded. 

We equip each pair (fc, a) with two labels, which are chosen from a finite set but which still 
locally distinguish between different pairs {k,a). The primary label, denoted by labcli(fc,a) S 
{0, 1, . . . where L is the maximal number of neighbours, is chosen in such a way that any 
two neighbours have a different label. The secondary label, denoted by labcl2(fc, a) G {1, . . . , M}, 
where M is the maximal number of children, is chosen in such a way that no two children of 
the same parent have the same label. These labels were introduced in [T^ with slightly different 
notation. 

Parametrized points and partial order. As described above, we now want to perform a 
perturbation of the original x'^ and < so as to obtain a parametrized family of similar objects, on 
which probabilistic statements can later be made. The parameter space will be 



({0,l,...,L}x{l,...,Af} 

with a typical point denoted by w = (a;/c)fegz, where ujk ~ {h, nik) <= {0, 1, . . . ,L} x {1, . . . , M}. 
The new dyadic points — z^{ujk) are defined by 



ix'^p'^^ if labeli(fc,Q!) = 4, and (fc -(- 1, /3) < (fc, a), and label2(fc -I- 1, /?) = m^. 



\x'^ if labeli(fc, a) 7^ 4, or ^(fc -t- 1, /3) < (fc, a) such that label2(fc -1- 1, /3) = mj,. 

A key feature of this definition is the following probabilistic statement, when D, is equipped with 
the natural probability measure P;^, which makes all coordinates wj, — {i^, nik) are independent of 
each other and uniformly distributed over the finite set {0, 1, . . . , L} x {1, . . . , M} (in other words, 
ik is uniformly distributed over {0, 1, . . . ,L}, and over {!,..., M}, independently). 

2.3. Lemma. Let (k < (fc,a) be fixed. Then 

In other words, every old point on the level k + \ has a positive (and bounded from below) 
probability of being a new point on the level k. 

Proof. A sufficient condition that = x^^^ is that 1^ ~ labcli(fc, a) (which has probability l/(L-f- 
1)) and that ruk ~ label2(fc -f- 1,/3) (which has probability 1/M). Multiplying the probabilities of 
these independent events gives the claim. □ 

The point of the next result is that the new points behave qualitatively like the reference points, 
only with slightly weaker constants. 



2.4. Lemma. The new points satisfy 

1 



d{zl,zl)>—6\ uAnd{x,zl) KAAlS"". 

ZAq a 



Proof. Consider the second bound first. Note that, in either of the two possibilities for the new 

mmd{x,z^) < mm[Aod{x,xl) + Aod{xl, z'^)] < 4Al6''. 



point, we have z^ = x'i'^^ for some (fc + 1,(3) < {k,a); in particular, d{x'^,z^) < 2AoS'^ , and thus 



Let us then estimate (i(z^,z^) from below for a ^ (3. If {k,a) and (fc,/3) are not neighbours, 
then by definition this distance is at least {2Aq)^^5'^. So suppose that these pairs are neighbours. 
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(2.5) + <^ (fc,a) 



Then they have different primary labels, and hence at least one of the new dyadic points, say 
0^, must satisfy = x^. On the other hand, if {k + 1,??) < {k,f3), then (fc + 1,77) ^ {k,a), 
and we thus know that , a;^) > {2Ao)~^S'' . But will be one of these points a;J^+^; thus 

d{z^,z^) > {2Aa)^^d'', as we claimed. □ 

The new partial order <i^, uj = {ujk)k&i is set up as follows. We declare that 

^ {d{xf\zl)<\A-H\ or 

\{k + l,(3)< ik,a) and ^7:d(x^+\z^) < ^A-^SK 

In other words, to find the new parent of (fc + for the new partial order <i^, we first check 
whether the reference point x^+i is close (within distance jAq'^S'') to some new dyadic point z^. 
If yes, then the corresponding (fc,a) will be the new parent of (fc + l,/3). If no such close point 
exists, then we simply use the original partial order < to decide the parent of (fc + 1, 

Properties of the new points and order. 

2.6. Lemma. For any given fc, the truth or falsity of the relation (fc + 1,/3) (fc, a) depends 
only on the component ujk oftu. 

Proof. In the defining conditions on the right of (|2.5p . the only dependence on to is via the new 
dyadic points z^, z^, and they depend only on Uk by definition. □ 

This explicit definition of <^ in terms of the geometric configuration of the points and the 
original partial order < is a novelty of our construction, where |19[ I21| required a condition similar 
to (|2.2p . which only specifies the relation up to certain degrees of freedom. For us, a condition 
analogous to (|2.2|1 is a consequence of the definition: 

2.7. Lemma. 

d(z^+\z^)<iV<5'= (fc + l,/3)<. (fc,a) =^ diz';+\z'^)<5AldK 

Proof Let first d(z^+\z^) < ^A-^S''. Then 

d(4+\z^) < A,d{x';+\z^+') + Aod{z^+\z^) < 2Al6''+' + ^A.^d'' < \A,'d'', 

and hence {k + 1, (3) <^ {k, a) by definition. 

Let then (fc + 1,^) <^ (fc,a). If d(x^+\z^) < jA-^S'', then 

d(^^+\ 4) < Aod(z^-+\ 4+1) + ^orf(4+\ 2^) < 2^25^+1 + < 5^3^^ 

Otherwise, we have (fc + 1, /3) < (fc, a), and hence 

diz';+\zt) < Agd(z^+\ 4+1) + A^d(4+\ 4) + Aodix^z^J 

< 2AlS''+^ + 2AlS'' + 2Al6'' < bA^S''. □ 

We can iterate this as follows: 

2.8. Lemma. For all £> k, 

d(z^,z^)<iV<5'= =^ (^,/3)<. (fc,a) =^ d{zf,,z'^)<6AtSK 

Proof. The second implication follows from the second implication of Lemma 12.71 with the triangle 
inequality: 

1 - /In/) ^ 

For the first implication. \i I = k, then the closeness of the points implies that 13 — a. li t > k, 
consider 7 such that {£,(3) <uj (k + 1,7). By what we just proved, d{zjj, z!^'^^) < GAqJ'^+i, and 
hence 

diz^+\zt) < Ao{lA^H'^ + 6At6'^+') < U^'S'^; 
thus {i,f3) <^ (k + 1,7) <aj {k,a), where the last step follows from Lemma [^771 □ 



(€, 13) (fc, a) =^ diz^zl) < bAl5'^+^ ■ A^o+i < < QAt& 
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Dyadic cubes. With the auxihary objects at hand, the dyadic cubes are easy to define. As in 
|19j . we introduce three families of these cubes — the preliminary, the closed, and the open: 

Qiiiv) :={2^H :(£,/?) <^ (fc,a)}, 

g^(L.) :=gSH, Q^(w) := interior Q^(l.). 

Note that Q^(aj), and hence (5a('^) a-nd Q^{ll!), only depends on u!e for £ > k. 

The rest of the section is concerned with the properties of these cubes. We first deal with 
several properties valid for an arbitrary fixed choice of the parameter a; S 51, and finally present 
a probabilistic statement concerning a random choice of lo under the natural product probability 
measure on the set f2. 

2.9. Lemma. 

Proof. Let x £ Then x is a limit of some points with {£,(3) <^ {k,a), and we may 

choose a subsequence with {£,P) <cj {k + 1,7) <uj {k,a) for some 7. Then 

d{z^^,x) < Aod{zlz';+^) + Aldiz';+\zj,) + A^dizj,, x), 

where the sum of the first two terms, by Lemma [2T7l and ([2]) is < ^AqS'^ + GAqS'^'^^ < 6Aq6'', and 
the last term becomes arbitrarily small. □ 

2.10. Lemma. We have the following covering properties for each fixed k G Z; 

Proof. Every x X is the limit of points z™^^_^-^ where m — > 00. When m > k, we have (to, 9{m)) < 
{k, a) for some a, and hence IJ^ <3q(w) is dense in X. On the other hand, the uniform separation 
of the centres z^, the uniformly bounded diameter of the sets Q'^{uj), and the geometric doubling 
property imply that this union is locally finite. Hence the closure of the union is the union of the 
closures. 

Similarly, it is immediate that Q'^{(^) ~ U0 (fc+i ^)< (k a) Q'p^^i'^)^ and the second claim follows 
by taking the closures, observing that the union is finite. □ 

2.11. Lemma. The closed and open cubes of the same generation are disjoint: 
In fact, 

QU^r = U QU^y 

Proof. We first check the weaker statement that 

(9^(a;)nQ^(w) = 

Indeed, for contradiction, let x G Q'^i'^) ^ Qpi^)- Thus x — zi^ with (£,7) <uj {k,P), and also 
X = lim,„_j,oo Zg^'^^j^ with {m,9{m)) < {k,a). For m large enough, we deduce that d{z^^^y z^) < 
and thus, by Lemma [2.81 {m,9{m)) 7) <uj a contradiction. 

To prove the actual first claim, again by contradiction, let x G Q^{oj) n Q^^{uj). Hence x = 
liuVm^ooXm with x„i G Since Q^p(uj) is open by definition and x„i — > x, it follows that 

Xm G Q^(i^) ^ for large enough to. So in fact Xm & Qq(w) n Qp{uj), but this is empty by 

the first part of the proof. The claim follows. 

For the second claim, it is immediate from the first claim that Q^i^) ^ r\/3^aQ'^(^y ~ 

( U/35^a Q'^i'-^)j ^a- On the other hand, by local finiteness, the union defining {^aY closed, 

and therefore itself, open. Since X = Q^iuj) U (ff^y, we have ff^ C and since Q^{uj) 

is the largest open set with this property, we get C Q'^{u}). □ 
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2.12. Lemma. 



Proof. It suffices to prove that B{z^, ^A^^S'^) is disjoint from each Qp{uj) with /3 ^ a. For 
contradiction, let x G B{z'^, ^A^^d'^) n Q^pito). Thus x ~ limm^oo ^ejm) for some {ni,9{m)) <^ 
{k,(3). Then 

d{zTi.n),4) < Aod(z,7„),x) + Aod(x,z^) < ^^A.H'^ 
for large enough m, since the second term is strictly smaller than this bound, and the first term 
tends to zero as m — > oo. But then Lemma 12.81 says that {m,6{m)) {k,a), a contradiction 
with {m,0{m)) <^ {k,P). □ 

The following theorem summarizes the above properties of the dyadic cubes for a fixed param- 
eter w, and supplements the key statement about their probabilistic behaviour under the random 
choice of e f2. 

2.13. Theorem. For any fixed uj € U := ({0,1,..., L} x {1,...,M})^, the cubes satisfy the 
following relations: the covering properties 

a I3:(k+l,l3)<^{k,a) 

the mutual disjointness property 

g^(a;)nO^(w) =0 (a^/3); 

and the comparability with balls: 

B{zl \A-H'^) C Qlioj) C Ql{^) C B{ziMl5^). 

Moreover, when J7 is equipped with the natural probability measure we have for some rj € (0, 1] 
the small boundary layer property: 

(2.14) F^[x e \Jd,Qi{uj)) < Ce" (a,Q^(c.) {y e Q^(c.) : %,^Q^(w)) < eJ^) ; 

a 

and in particular the negligible boundary property: 

k,a 

Henceforth, r] will always designate the fixed positive constant provided by this 
theorem. It will reappear as the Holder-regularity index of our splines and wavelets. 

Proof. It only remains to check the probabilistic statements. 

Let X £ X , k E 'Z, and £ > be fixed. For every i = k, k + 1, . . ., there is some 7 so that 
d{x,xi^'^^) < 2AoS^^^. Now suppose that for some i, the point x^"^^ is chosen as a new dyadic 
point z^, and recall that B{zjj, ^Aq^S^) C Q^^iu}). Thus, 

d(x,^Q^(c.)) > ^^d{zl^Q'p{oj))-d{x,/p) > (i^-e - 2Ao6)6' > l^A^'S' > eS^ 

provided that i < k + log(7j4ge)/ log(5. In particular, we have x £ Qp{i^) C Q'^{u!) if <uj 
(fc, a), while d{x^ ^Q^{uj)) > d{x, ^(5^(i^)) > e^^^. Thus x cannot be in IJ^ d^Q^{uj) in this case. 

In other words, in order that x belongs to the union [Ja'^^Qai'^)j it is necessary that none 
of the points a;^+^, where k < £ < k + log(7AQe)/log(5, is chosen as a new dyadic point. But, 
by Lemma 12.31 every xi^~^^ has a probability at least t := j^jj^^^i-^ of being chosen. So the 
probability that a;^+^ is not chosen is at most 1 — r. Moreover, these events for different levels 
£ are independent of each other. Hence the probability that none of the 1^+^ is chosen, for 
k<£<k + \og{7Ale)/ \ogS, is at most 
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where r/ log(l — T)/\og5 > (since both 6,1 — t G (0, 1)) and C = (TAg)''. This is exactly as 
claimed. 

The negligible boundary property follows from the small boundary property as e — > 0. □ 

We conclude this section with the observation that in our construction, the original dyadic point 
may also be viewed as a 'centre' of the new dyadic cubes 0^(0;) for all w G fi: 

2.15. Lemma. B{x'^, ^A^^S'') C g^(cj) C i?(a;^, 8AgJ'=) for any uj eQ. 

Proof. By Lemma [2.91 we have for any a; e that 

dix, xi) < Aod{x, zi) + Aod{zi,xl) < GAlS'' + 2AlS'' < SA^J^ 

which shows that C _B(x^, 8Ag5'=). 

For the other inclusion, let d{x, x^) < ^A^^S''. Since X = IJ,, Qf-i'^^i^)^ we have a; € Q,')"'"^'^) 
S(a;J;+\8Ag(5''+i) for some ?/. But then 

dix';+\xl) < Aod{x';+\x) + Aodix,xi) < 8Al5''+' + i^^' < K"'^'- 

By another application of (|2.2p . this implies that {k + 1,1]) < (fc,a). 

We want to prove that in fact (fc + 1, 77) (fc, a), since then a; € Q^'^^{uj) C Q'^(lj). By (|2.5p . 
the only potential obstacle to this is that d{x^'^^ , z!^) < jAq'^S^ for some 7 ^ a. So suppose this 
is the case, and recall that z^^ = ^e^^ for some (fc + 1,0) < {k, 7). Then 

< Aod{z';,x'^+') + Aldix';+\x)+Ald{x,x':,) 

and hence (fc + 1,0) < {k, a) by (|2.2p . But also [k + 1,9) < {k, 7), hence 7 = a, a contradiction. 

We have shown that an arbitrary x € B{x'^ , ^Aq^S'') satisfies a; £ Qj^^^ (w) for some 77) 
{k, a); thus X e Q^{uj). □ 

3. Construction of splines 

The construction of splines on X, and the proof of their basic properties, is now an easy 
consequence of the preparations from the previous section. For every (fc, a), we define the spline 
function 

(3.1) slix):=F^(^x€Qi{Lu)). 

3.2. Theorem. The splines (|3.ip satisfy the following properties: bounded support 
(3-3) 15(^.^1^-3^^) (2;) < siix) < lB(x^,8A=5'')(a;); 

the interpolation and reproducing properties 

(3.4) s^(a:^)=<5„^, ^4(a:) = l, s^(:^) = ^P^^ • 4+^(x) 

a p 

where {p^p}p is a finitely nonzero set of nonnegative coefficients with "YlipP^p ~ 1/ '^^'^ Holder- 
continuity 

\sii^)-siiy)\<c{^)\ 

Regarding p.3p , remark that it is a bit unusual that a spline be a non-zero constant on part of 
its support. 



Proof. The relations p.3p are immediate from Lemma 12.151 and this implies in particular that 

si{xi) = 1. 

Since the boundaries have vanishing probability and X = IJ^ Q^{lu), it follows that 
(3.5) ^s^(x)=pJa:ex) =1. 
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Since these functions are nonnegative and sjj(x„) = 1, it must be that s^(x^) = for /3 7^ a, and 
hence in fact one has the interpolation property 

(3.6) 

From the identity 

/3:(fc+l,,3)<^(fe,a) 

(and using again the vanishing probabihty of the boundaries) we also have 



I3:(k+lj3)<^(k,a) 



^P„({(fc + i,/3) (fc,a)} n [x e Qf\Lo)] 



13 

= ^P^((fc + l,/3) {k,a))sf\x) Y^p^^ ■ s)+\x), 

13 f3 

where the key third step used the independence of the two events; namely, the event (fc + 1, /3) <^ 
(fc, a) depends only on oj^-, while the cube (5^^^(w) depends on for £ > fc + 1. The support 
properties of the splines readily imply that only boundedly many of the coefficients are nonzero 
for a given (fc,a), so that span{s^}Q C spanjsj^'''^}^. 

The Holder-continuity of the splines follows from the probabilistic smallness of the bound- 
ary regions, as expressed by p.l4p . Indeed, 



< 



< P. 

< c 



d{x,y )\^v 
gk 



□ 



4. Auxiliary results on quasi-metric spaces 

We include this intermediate section to collect some auxiliary material so as to streamline 
the subsequent analysis. This section is mostly concerned with difficulties of a quasi-metric in 
contrast to a metric. The results are mostly part of the folklore, but somewhat difficult to find 
in the literature in full generality, since additional assumptions on the space are usually imposed. 
(After completion of this work, we became aware of the recent preprint |31| where the density 
result below is also proved in full generality by a completely different argument.) 

For the beginning of this section, as above, we only assume that X is a geometrically doubling 
quasi-metric space. Throughout, let 77 > denote the fixed positive constant from the 
small boundary property of the dyadic cubes and the Holder regularity of the splines. 
4.1. Lemma. Let F C G C X be sets with 



diF,G''):= inf fi(a;, y) = A > 0. 
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Then there exists a Junction : X — > R with Ij^ < "yS < 1g '^'^c' 
Proof. Let k be the smallest integer so that 1%A15^ < A. We set 

Then f{x) = 1 for x £ since ■^aC-'') = 1' ^^^^ ^^^^ sum defining (p{x) contains all s^(a;) 
whose support intersects F. To prove that (p{z) = for z £ G"^, we need to check that no 
suppsj^ C B{x^,8Al5'') appearing in ip can intersect To this end, let y G B{x'^,8AlS'') Ci F 
(which exists by definition) and z G B{x'^,8Aq6'') C (whose assumed existence should lead to 
a contradiction). Then 

d{F, G^) < diy, z) < Aod{y, x^) + Aod{xi, z) < 2Ao ■ SA^S'' < A, 

a contradiction indeed. So we have If < V Id a-nd it remains to check the Holder-continuity. 
By geometric doubling, there are only boundedly many s^, whose support contains either x or y. 
Each of these satisfy the estimate 

\sii^)-siiy)\<c{^y<c{^)\ 
and hence so does their sum over boundedly many indices a. This completes the proof. □ 

4.2. Corollary. For every B{x, r), there exists a function ip : X R. with lB{x.r) l£ ^ l£ ^B{x,2Aor) 
and 

\^ix)-piy)\<c{^y. 
Proof. If y E B{x, r) and z G B{x, 2AorY, then 

d{y, z) > ^rf(z, a-) - d{x, y) > 2r - r ^ r; 

^0 

thus d{B{x, r), B{x, 2AQrY) > r, and the previous lemma applies. □ 

We formulate a quasi-metric version of a well-known covering lemma for metric spaces, cf. [181 
Theorem 1.2]; the quasi- metric extension is obtained mutatis mutandis. 

4.3. Lemma. Let he a family of balls B{x, r) in X , with hounded radii and contained in a 
hounded set. Then there exists a pairwise disjoint suhcollection {Bi\^^, whose concentric expan- 
sions hA^Bi cover all original B{x, r) G 

4.4. Lemma. X is a countable union of hounded open sets. 

Proof. Recall that the interior of a set E is the open set {E'^Y C E. It is easy to check that 

B{x.,r/Ao) C interior _B(a;, r); 

hence, with any base point xq G X, the bounded open sets interior i3(a;o, n), n G N, provide the 
required covering. □ 

So far everything has been based on the space geometry, i.e., the properties of the quasi-distance 
only and geometric doubling. We next add a measure /i into our considerations. 

4.5. Proposition. Let ^ he a non-trivial Borel measure on X , finite on bounded Borel sets. Let 
1 < p < oo. Then H6lder-iy continuous functions of hounded support are dense in L^{^), where rj 
is the Holder exponent of the splines. 
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Proof. By the density of simple functions due to general measure theory, it suffices to show that 
for every bounded Borel set E and every e > 0, there exists a boundedly supported Holder-ry- 
continuous function (p with \\1e — f\\p < 2e. By a general result concerning Borel measures [121 
Theorem 2.2.2 (ii)0, there is an open set G ^ E such that /Lt(G \ E) < e^. For every x G G, we 
choose a radius small enough so that B{x,AAlr.j;) C G. The balls B{x,5-^AQ^r^) cover G. 
and hence there are pairwise disjoint balls B{xi,5~^AQ^ri), := r^.., so that B{xi,rj) cover G. 
As we have B{x,r) C B[x,r) C B{x,2Aor) for any ball, with B being the topological closure of 
the baU B, 

OO OO 

G C y B{x,,r,) C y 2Aon) C G. 

i=l i=l 

Thus 

yu(G) = lim /i( I \B{xi,ri)] =: lim ju(G„). 

71— >-C30 \ ^-^ / n— >cjo 

Let us fix n so large that /i(G \ G„) < e^. We check that (i(G„, G"^) > 0. Indeed, if a; € G„, then 
a; € B{xi,ri) C _B(a;i, 2Aori) for some i = 1, . . . ,n, whereas B{xi,4:AQri) C G. Hence 

d(G„,G'=)> min d(B(.T,;, 2Aor,;), B(a;,, 4^2?-,)=) > min r, > 0. 

z— l,...,n ?'— l,...,n 

By Lemma |4.1[ there exists a Holder-7/-continuous with 1^^ ^ ^ ^ 1g- Hence ||1g — ^\\p ^ 
II 1g - IgJIp = G„)i/P < e, and finally ||li, - ^\\p < 2e. " □ 



5. L' THEORY AND MULTIRESOLUTION ANALYSIS 

We now return to the development of the spline theory in the presence of a nontrivial Borel 
measure fj, on {X,d). From now on, we assume that (X, d, /i) is a space of homogeneous type in 
the most general sense, namely that in addition to having a quasi-metric d we only assume the 
doubling condition which reads: for all x E X and r > 

< n{B{x, 2r)) < Cf, ■ ^i{B{x, r)) < oo. 

This inequality makes sense if balls are Borel sets. But this may not be the case. However, 
their topological closures are (because they are closed sets by definition) and satisfy B{x,r) C 
B{x,r) C B{x,2Aor). In that case, up to changing constants throughout, we may change balls 
to their closures in the above definition of the doubling condition, thus avoiding to resort to the 
outer measure associated to fi. To simplify matters though, we assume that balls are Borel sets 
and leave to the reader the modifications if not. 

The doubling condition on (X,d,^) implies that {X,d) is geometrically doubling (see [5]), 
so that the spline functions, built independently of measure considerations, are at our disposal. 
We now show that they provide a multiresolution analysis of i^(^), essentially in the sense of 
Sweldens [32! Definition 3.1], who considered the case of general measures on M". This consist 
of all properties of a classical multiresolution analysis of Meyer [301 Definition 2.1], to the extent 
that this definition is meaningful in a quasi-metric space context: the classical postulates dealing 
with translations and dilations, specific to the Euclidean space and the Lebesgue measure, are 
now meaningless. 

5.1. Theorem. Let Vk he. the closed linear span o/js^jja in L'^{fj.). Then Vk C Vk+i, and 

{0}, if X is unbounded, 

Vkg ~ {constants}, if X is bounded, 



\JVk = L^^,), f]Vk = 



This result is stated in a metric space with an outer measure but the argument is valid for any topological space 
for a Borel measure. Moreover, it is assumed that E can be covered by countably many open sets with bounded 
measure which is granted from Lemma 14.41 and our assumptions on /j,. 
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where /cq is some integer. Moreover, the functions s^j/y^T*^ form a Riesz basis of V^: for all 
sequences of numbers Aq, we have the two-sided estimate 



1/2 



with ^* ^l{B{x'^,S'')) 



Proof. The nesting property Vk C 14+ 1 is immediate from the reproducing properties of the 
sphnes. 

The Riesz basis property. We use properties (j3.3|) and p.4p of the sphnes. At any point x, the 
values Sa(a;) are nonnegative and sum up to 1 (when a is the summation variable and k is fixed). 
Hence 

CK a a 

Thus 

f I ^ A„4(x)|' dM(x) < ^ |A„|V(B(xL SA^^'^)) < |A„|V(S(2:;;, <5''0) 

a a a 

by the doubling property. 

On the other hand, we also know that Sp{x) = 1 for a:; € B{x'p, g^i^'^) and all other s^, a ^ (3, 
vanish on this set. Hence 



thus 



and therefore 



EA"^^(^) >Ei^"i'ii 



I |^A„4(a;)|'dM.T)>El^"lVB(-^o,8ij'5'=))>l^|A„|V(i3(x^,'5^^ 
again by the doubling property. 

The union of the Vfc. If / is a Holder-ry-continuous function with bounded support, the sum 
fk{x) = J2a f{^a)^ai^) defines an element of Vk and using p.4p we have 

/(x)-/.(x)=5](/(x)-/(x^)).^(x) 



so that 



||/-M|oo<sup sup \f{x)-f{x',)\<C6'^. 

a xeB{x'^,8AlS'') 



Convergence in (/i) follows from this and bounded support for f — fk- 

The intersection of the Vk- In case X is bounded, it is clear from definition that the sets {x^ja 
reduce to one point when k gets small. In that case, s^{x) = 1 for all x e X by (13. 4[) . So we 
choose fco as the largest integer such that this property holds. 

If X is unbounded, then ^i{X) = oo as /i is doubling. (See remark after Lemma l8. II for a proof 
of this known fact.) If / G HfeGZ ^fe' ^^'^'^ II/II2 ^ J2a for all k. Consider a fixed x € X. 

Then f{x) = J^a fi^a)^a{^)i where the sum is over the boundedly (with respect to k) many a 
such that a;^ € B{x,8A^,S'^). Thus |/(a.)| < E J/(4)l < E JI/lb/^M^ < ll/IUM^, <5'=). This 
tends to zero as fc — ?• — 00, and hence f{x) = 0. Since the point was arbitrary, we have / = 0. □ 
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6. BlORTHOGONAL AND ORTHOGONAL SPLINE SYSTEMS 

In this section, we use a classical algorithm (cf. [3D1 Sec. 2.3]) to construct two further bases of 
the space Vk spanned by the splines {Sq}q. We use the abbreviations 

for these frequently appearing volumes. 

6.1. Theorem. There exist a system of hiorthogonal splines {.5^}q. in Vk with 
(6.2) {si,tp)LH^)=5o.,0. 

as well as an orthonormal basis {^^ja ofVk, which satisfy the following estimates, where d{x, y) < 



l^i ■ + Va'S • |</'^(x)| < Cexp ( - ^5-^d{xi,xf ) 



with s = (1 + 1052^0)^"'^ or s = 1 if d is a Lipschitz- continuous quasi- distance. 

6.3. Remark. From the decay estimates and doubling, it readily follows that g L^{iJ.). Summing 
the biorthogonality relation (|6.2p over all a and recalling that J2a ■^q(-'^) = '^^ deduce that 



J dti = l. 



We begin the proof of Theorem 16. II Fix /c e Z. By abuse, we identify = {a;Jj}Q, and the set 
of indices a corresponding to points a;^. The previous section says that there is a linear, bounded, 
injective map Uk - L^if^) with closed range, defined by 

for A = {Aajjjg^fc. Here, £'^{S''') is equipped with counting measure. Let Vk denote its range. 
By the properties of the splines, the sum defining each element of Vk converges locally uniformly 
and defines a Holder-Ty-continuous function. The inverse of Uk can be computed using (j3.6p by 
{Uk'fh = fixi)^. 

Let (5^ be the canonical orthonormal basis element in £^(^'''). By construction UkS'^ ~ ^° 



Since Uk is an isomorphism, this means that the splines sjj form an unconditional basis of Vk. To 
find a biorthogonal system in Vk, that is a system such that 

we observe that if / = f/fcA, /' — Uk^ then 

(/, ^ (^4-A, A')£2(^fc) 
with Mk being the infinite matrix with entries 

Mk{a,l3) = . 

Another reformulation of the isomorphism property of Uk is that Mk is bounded and invertible 
on ^^(j£"'°). It is also positive self-adjoint. So the biorthogonal system is uniquely defined by 

M^s^ = UkM^'6l 
If one wants an orthonormal basis of Vk C L^(/^), one defines instead 

= UkM7^/^5l 
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In other words, 

We now establish decay estimate of the coefficients Afjr^(a, /3) and Afj, ^^^(a,/?). 

6.4. Proposition. There exist constants 7 > and C < 00, independent of k such that the entries 

■I 1/2 

of M^^ and ' have upper bounds 

Cexp(-7(r'^d(a,/3))^), 

with s = (f + log2 Ai))^^ or s = 1 if d is Lipschitz- continuous. 

If we introduce the induced normahzed quasi-distance on ^ x dk{a,f3) := — fr^, we 

have to prove uniform estimates on the entries of M^T^ or Mf, in terms of dk- Note that for 
this quasi-distance is a band-matrix, more precisely Mk{a, /3) = if dk{oi, /3) > 16 Aq by (|3.3p . 
hence it has exponential decay. 

6.5. Lemma. Let 'E. be a 1-separated set in a quasi-metric space (X, d) with quasi-triangle constant 
Aq having the geometric doubling property with constant N . Then for all £ > 0, there exists 
c{eT Aq^N) < 00 such that 

sup exp (^sd{a, E)/Ao) E '^'^P ( ~ £d{a, /?)) < c(e, Aq, N). 

Moreover, for any < c' < c/Aq 

E 6xp ( — cd{a, P)) exp ( — cd(/3, 7)) < C exp ( — c'd{a, 7)) 
/3eH 

where C does not depend on a, 7 G X. 

Proof. Fix a point a G X. Pick any point G S. We have 

E exp ( - ed(a, /?)) < exp ( - £d(a, /3')Mo) E ^xp ( - erf(/3', /3)Mo) . 

/3GH /3GH 

The number of points in S at quasi-distance at most 2-' from j3' is bounded by C-'+^ from some C 
depending only on and N but not /?'. Thus 

Eexp(-£d(/3',/3)/Ao) <C-hE E exp(-£d(/3',/3)/Ao) <E^-''+'-exp(-£2VAo). 

/3GS: j>1 /3:d(^i',/3)~2J i>0 

Taking the infimum over all /3' € 'E, proves the desired inequality. 

The result for the matrix product coefficients follows analogously and we skip details. □ 

6.6. Lemma. Let E be a l-separated set in a quasi-metric space {X, d) with quasi-triangle constant 
Aq having the geometric doubling property with constant N. Consider a matrix M = (M(a,/3)) 
indexed by E x E such that there exists c > for which 

C= sup exp (c(i(a,/3))|Af(a,/3)| < 00. 

Then M is bounded on ^^(S). If M is invertible, then there exists c' > such that 

sup exp(c'(d(a,^))^)|A//-^(a,^)| < 00 

(",/3) 

with s = (1 -|-log2 Aq)^^ or s = 1 if d is a Lipschitz- continuous quasi-distance. If, in addition, M 
is positive self-adjoint, then the same conclusion holds for M~^l'^ . 
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6.7. Remark, (i) Note that the exponent s = 1 for a usual distance is recovered as a special case 
in two ways, either by setting = 1 or using the Lipschitz-continuity. 

(ii) The exponent s is in general optimal in this result. Namely, consider the band matrix 
M indexed by Z with M{i,i) = 1 and M{i,i + 1) = -A e (-1,0) for aU i G Z, and all other 
entries equal to zero. Then M~^{i, j) = A'~^ if i > * and zero otherwise. Now equip Z with the 
quasi-distance d(i,j) = \i — where r > 1. Its quasi-triangle constant is Aq = 2'"^^, and hence 
the exponent given by the lemma is s = (1 + log2 Aq)^^ = (1 + r — 1)^^ = 1/r. This gives precisely 
the correct decay of the matrix with respect to d, as 

A'"-' = cxp(-c • \i - j\) = exp(-c • d{i,jy^''), c = log A~"^. 

Proof. If d is a genuine distance or a Lipschitz continuous quasi-distance, then this follows from 
Theorem 5 in [53] and the remark that follows it, which extends an earlier result in jTD] for 
band- limited matrices (but we shall need the full result here), with s = 1. However, as said not 
all quasi-distances are Lipschitz- or even Holder-continuous and we provide an argument in full 
generality, which also recovers the mentioned special cases. 

We begin with the following observation. For n> 1, let k„ be the best constant in the inequality 

d(ao, a„) < K„(d(ao, ai) -I- d{ai,a2) + . . . + (i(a„_i, «„)) 

for every chain (ao,ai, . . . ,«„) of n + 1 elements (not necessarily distinct) of S. It is clear that 
(k„) is non-decreasing and ki = 1, K2 < ^o- Moreover, using d{ao, a„i+n) < A(){d{aQ,am) + 
d(am, «„_(_„)), it follows that Hm+n l£ ^0' nrax{KTO, K„}. Thus Aqk^ and therefore, < 

A^. We conclude that k„ < Aq"*"'"*^^" = ^0"-'°*^"^°- Note also that if d is i-Lipshitz, then 
d{ao, a„) < d(ao, an-i) + Ld{an-i, Q;„), which by iteration gives Kn < L for all n. 

Now assume that M is positive self-adjoint and invertible. In this case, one can write M = 
h{I — A) with h = {\\M\\ + |jM~^||)/2 a positive real number and A a matrix with norm r = 
{\\M\\ - ||M-i||)/(||M|| + |lM-i||) < 1. Moreover, the coefficients of A have the same decay 
as those of M. Without loss of generality, we normalize h = 1. Develop {I — A)~^ in the 
Neumann series ^A^ and estimate the coefficients A^{a,(3), n > 1, a ^ /3, in two ways. First 
|A"(a,/3)| < r". Second, we have 

\A"{a,P)\< C"cxp{-c{d{a,ai)+d{a,,a2) + ...+d{an-i,l3))) 

(ai,...,Q„_i)eH"-i 

< C" exp ( - -^d(a, ^)) ^ exp(- ^(d(ai,a2) + ... + d(a„-i,/3))) 

(qi,...,Qti-i)G^"--^ 

<C"cxp(--^d(a,/3)), 

where we applied Lemma 16.51 rt — 1 times with e = c/2, and C = C ■ c(e, Aq, N), in the notation of 
that lemma. As k„ is non decreasing, we have for any integer no using the second estimate for 
< n < Hq and the first for n > Hq, 

\M-^{a,f3)\ < (no + l)C'"«exp(---^d(a,/3)) + r"°+i(l - r)-i 

and {no + 1)C'"° < for some large constant D > 0. Choosing uq as the first integer such that 
the ffist term dominates, we see that d{a, j3) ~ no • k„o ^ "■0^'', where s = 1/(1 -l-logj Aq) (or 5 = 1 
if d is Lipschitz-continuous, recalling that k„q < L in this case). Hence, for some constant c' > 0, 

\M~^{a,l3)\ < r"° < exp ( - c'd(a, ^)" ) . 

For A/^^/^, we use the power series (/ — A)~^^^ = ^c„A". As < c„ < 71^/^, the argument 
is the same. Finally, if M is not positive self-adjoint, then we use M^^ = M* {M M*)^^ and the 
remark that on a 1-separated set d> d" . □ 

Proof of Proposition \6.4] and Theorem \6.1\ Lemma 16.61 gives us the desired decay for M^^ and 
Mj. ' , as claimed in Proposition l6.4[ with uniform control of the constant with respect to k. The 
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decay and regularity of the sjj and 0jj , as asserted by Theorem 16.11 then follow from the support 
and regularity of the and the decay of the matrix coefficients. We skip details. □ 

6.8. Remark. Let < s < 1. If d is a quasi-distance with quasi-triangle constant Aq, then is 
a quasi-distance with quasi-triangle constant Aq. If d is geometrically doubling with constant N, 
then so is with constant at most A^Ti/sl^ where [ ] is the 'rounding-up' to the next integer. If 
/i is a doubling measure with respect to the balls for d with constant c^, then it is with respect 
to the balls for with constant at most cj}^^'^ . Thus both previous lemmata apply to d^ with 
constants that depend only on the original constants and s. We shall use this remark later. 



7. Spline wavelets 

We are now prepared for the construction of an orthonormal basis of consisting of 

wavelets ip'^ with similar decay and regularity properties as with the spline systems. We follow an 
algorithm from Meyer |29| . 

Fix k E Z, k > kg, where feg is as in Theorem 15. II Recall the operator 



is an isomorphism from Vk+i onto £'^{^'''^^). Denote by Yk the inverse image of the subspace of 
sequences vanishing on Jl^'^; this subspace is naturally identified with £'^{'3^'^), where 
we recall the notation := JT'^+i \ Clearly, 14 ® Yfc = Vk+i topologically. Consider the 

orthogonal (in L^(/z)) complement Wk of 14 in Vk+i- Then the restriction to Yk of the orthogonal 
projection Qk onto Wk is an isomophism onto Wk- Now, identifying with '3^'^ the set of indices 
/3 corresponding to x'^'^^ £ S(^^^^ \ = '3^^, the collection forms an unconditional 

basis of Yk and its image under Qk is an unconditional basis of Wk- A representation of Qkf when 

/ e Vk+i is 

(7.1) Qkf^f- J2 if^Si)LH^.A^f- E if^'')LH^.)sl 

because the sum is the orthogonal projection of / onto Vk- Hence, the pre- wavelets 

r,fc/ \ ^ k+i 



rp{x) ■.= Qksl+\x) 



have the L°°-normalized exponential decay Cexp ( — "f{6 ^d{y^,x)y ), where — ^^+1 



(Remark that they are normalized as the splines s*j^^.) Finally, one can orthonornialize them in 
L^(yu) following the procedure of Section IH] applied to the positive self-adjoint matrix 



M(a,/3) 



indexed by x that is. we define 



Note that with the notation = we have = ^(S(?/^, (5*^+^)) - fi{B{y^p, S'')). The point- 

set W'^ is a 1-separated set for dj,{x,y) := {d~^~^d{x,y)Y and the matrix M has the exponential 
decay, as required in Lemma with respect to d^. By Lemma and Remark 16.81 M~^/^ has 
decay 

< exp ( - y4(4, 4)' ) = cxp ( - 7(rt(x;;, x^))- ), 

where s = {1 + log^ A'^)-^ = (1 + (1 + log2 ^o)"Mog2 Ao)-i - (1 + loga Ao)(l + 2 log2 ^o)"^ or 
s = 1 if d is Lipschitz-continuous. 
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This yields an ortho normal basis tl;^{x), G of having the normalized exponential 
decay 

\ip^pix)\ < =cxp(-7((5-''d(4,x))°), a = sS. 

Gathering the construction for all k (and adding constants if X is bounded), one obtains an 
orthonormal basis of made of spline wavelets. 

7.2. Theorem. Let {X,d,ii) be any space of homogeneous type with quasi-triangle constant Aq, 
and a := (1 + 2 log2 Aq)^"'^ or a :~ 1 if d is Lipschitz-continuous. There exists an orthonormal 
basis ?/;^, fc G Z (and k > if X is bounded), y^ G '3^'^, of L^{fi) (or the orthogonal space to 
constants if X is bounded) having exponential decay 

< ^ exp(-7(^-"%g,^)r), 

Holder-regularity 

|V'^(x)-^^(y)|<^^^=(%^)''exp(-7(rt(y^,x)r), d{x,y)<S\ 
^^i{B{y%5'^))^ ^ ' 

and vanishing mean 

[ ip''p{x)dn{x) =0, keZ,k>ka,yle'W^. 

Proof. It remains to see vanishing mean and regularity. 

Using / s^dji ^ 1 and J2a = ^- deduce from (|7.ip for / G that 

I Qkfdfi = I /d^-X^(/,s^) I S^d/i = I /dA^-E / Md^ = / /d/i-/ /dM = 0. 

Since the pre- wavelets ?/;^ lie in the range of Qk by definition, we have J ip^dii = 0, and the same 
result for the wavelets follows from the convergent series representation of in terms of the 

As for the regularity, recall the smoothness of the splines: |s^(x) — s'^{y)\ < C{d{x,y)/ 5^)^ . 
This implies for the pre-wavelets the estimate 

IV-^Cx) - my)\ < - + Y.Ks'+'rsDWslix) - sliy)\ 

a 

<c{^^)\siy^^^^s^)i^)+Y^Cc^v{-^idiy',,xi)/S''r^ 

a 

<c{^y e^p{-^{d{x,y',)/S'r). 

For the wavelets, finally, we have 

„ exp(-7(%^,y^)M'^+i)'') 
^ ^ti{B{y^^,5'^)) 

< 



7 y^s^^ 



Note that the value of 7 > changes from line to line in these computations. Also we used that 
a < s and a variant of Lemma 16.51 for the quasi-distance d{x, y)l5^ on X. □ 
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7.3. Remark. The construction and also the next sections suggest that the label k + 1 would be 
more appropriate than k for Wk, the wavelets and their scale (5*^, because its keeps closer to 
the definition of the point sets a subset of which will take an important role. We have 

kept the wavelet community notation as in |30| . 



8. Technical estimates related to vanishing annuli 

We break the development of the wavelet theory with this technical section, which will provide 
us with useful estimates to streamline the subsequent presentation. A basic difficulty related to 
general spaces of homogeneous type, as opposed to those with the reverse doubling property, is 
the possible existence of arbitrarily large empty annuli B{x,R) \ B{x,r) = 0. This leads to a 
certain dichotomy: locally, we either have the reverse doubling estimate, or the vanishing of a 
certain annulus, both of which provide certain control, which we need to exploit in different ways. 
This is quantified in the following. (The next result is well-known, cf. |141 Remark 1.2], but we 
include it here for completeness, since we need it to derive some consequences which appear to be 
new.) 

8.1. Lemma. For every .t € X and R > r > 0, at least one of the following alternatives holds: 

fi{B{x,R))> {l + e)n{B{x,r)) or B{x, -^R)\ B{x,2Anr) = 0, 

ZAq 

where s 1/C^(3Aq), and C^(t) is the smallest constant such that n{tB) < C^{t)ii(B) for all 
balls B ex. 

In particular, as i? — > oo, we observe the following (again well-known fact): If there is a ball 
B{x,r) such that ij{X) < (1 + Cf,{3Al)-'^)^i{B{x,r)), then X = B{x,2Aor). If niX) < oo, such 
a ball always exists, and hence diam(X) < oo. 

Proof. Suppose that the annulus is nonempty, and let y € B{x, jx^R) \ B{x,2AQr). Let p := 
^j^d{x,y) > r. We claim that 

Biy,p)CBix,R)\B{x,r). 

Indeed, if z G B{y, p), then 

d{z, x) < Aod{z, y) + Aod{y, x) < Aop + Aod{x, y) < R, 

while 

d{z, x) > -r-d{y, x) - d{y, z) > -r-d{x, y) - p = 77-5-^(2;, y) > r. 

Ao Aq ZAq 

We also claim that 

B{x,r) C B{y,3Alp). 

Indeed, if w G B[x, r), then 

d{w, y) < Aod{w, x) + Aod{x, y) < A^r + 2Alp < 3Agp. 
Now it follows that 

p{B{x, R)) - ^i{B{x, r)) > M(B(y, p)) > c^^^j^^f i^iy, SA^p)) > ^-l^p{Bix, r)). □ 

The following lemma relates the mentioned dichotomy to the distribution of the dyadic point 
sets ^'=: 

8.2. Lemma. For every x X and r > 0, there exists a decreasing sequence, finite or infinite, of 
integers {kj}j^Q such that r < S'"' < 5^'^ < . . . such that 

V{x,S'')>{l + eyV{x,r) and d{x,'3^^) + 5^ > 5^^+^ if k-j > k > kj+i, 
where we intepret kj+i :— —00 if J < 00. 
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Proof. Let k{0) be the largest integer with ^'^'"^ > r, and let k{j + 1) be the largest integer with 
V{x, J'^O'+i)) > (1 + s)V{x, J'^^^'). (Note that the sequence terminates if and only if ^J,{X) < oo.) 
Thus < {l+e)V{x,6''^^'^), and hence 

2Ao 

For fc > fc(j' + 1) + 2, the ball on the left contains at least one element of For k < k{j) — 1, the 
ball on the right contains at most one element of J^T'^. Since the balls are equal, for k{j + 1) + 2 < 
k < k{j) — 1, the ball contains exactly one element of So if k and k + 1 are both in this range, 
i.e., if k{j + 1) + 2 < fc < k{j) — 2, then the intersections of the ball with 3^'^ and ^'^'^^ coincide; 
hence there is no point of ^T*^ = JT'^+i \ in the ball. Thus d{x, '3^^) > ■ J'^O+i). On the 
other hand, if fc S {k{j + 1) - 1, k{j + l),k{j + 1) + 1}, then clearly S'' > 6 ■ 6''^^+^\ and hence 

d{x, '3/^) + 6^ > T-r5''^^+^'> if k{j + l)-l<k< k{j) - 2. 

Also, for k in the same range, we have 

V{x,5^) > V{x,5^^^^) > {l+eyV{x,5''^°y) > {l+eYVix,r). 
So the claim follows by relabeling kj+i := k{j) — 2 for j G N and fco k{0). □ 

Sums of the following type appear in connection with the wavelets: 
8.3. Lemma. For all x ^ X and r^v,a > 0, we have 

V{x, 5^)-" cxp ( - i{6-^d{x, 3/^)^ ) < V{x, r)-r 

k:SI'>r 

Proof. Let kj be the sequence as provided by Lemma [8T2l Then 
J2 V{x,6'')-''exp{--f{S-''dix,?V'')y) 

J 

E v{x,d')-''cxp{-j{5-''d{x,3^'')r) 

j=0 k:kj>k>kj-ii 
J 

<E E (l+e)"'""t^(a;,0""exp(-7(5'^(^-^+i-'=') 

J— k:kj>k>kj-\-i 

J oo 

< V{x, r)-'' ( ^(1 + s)-^") ( cxp(-7<5-™")) < V{x, r)-^ □ 

J— m—l 

9. Technical estimates involving the wavelets 

Actually, the estimates here are valid for any family of functions -0^ which satisfy the same 
size and regularity estimates as the wavelets. This includes the condition that ifi^ be concentrated 
around the point S '3^^ , and the structure of the point sets '3^^ is important for some of the 
following estimates. 

9.1. Lemma. Let a as in Theorem \7.2[ For a fixed fc G Z, 

E \^'i^)€iy)\ < ySw) ( - li^^'dix, iV^))^ ) exp ( - i{5-^d{x, y)Y), 

and, for d{x,x') < ^j^d{x,y), 

aGAfc 

C . r fdix,x')\v 



< 



V{x,5'') I ' V d''^ 



{l' {^^^Y} ( - ^i^~'d(x, 3^^)T ) exp ( - i{5-H{x, y)r). 
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Note that if. instead, we have a family of functions if^ corresponding to the points G 
instead of € '3^^ , we get exactly the same estimate but with d{x, ^*') in place of d{x, '3^^) in 
the result. Since ,'^'^) < 2Ao6'', the first exponential factor is roughly 1, and may be dropped. 

Proof. By the doubling condition and the quasi-triangle inequality, we estimate 

< ^4tt ( - ^(^~'d{x, vDT ) exp ( - i{5-U{v, y^T ) 

- V^W) ""^^ ^ ~ 2^'))'' ) exp ( - l{5-''d{x, y)r ) , 

and the sum over a G ^V'^ of the second factor is dominated by exp ( — ^{5^^d{x, '3^^)Y )■ 
If 5^ > d{x,x'), then 



< 



C fd{x,x')\v 

(^%^)" <5^P ( - li^^'di^x, yi)r ) exp ( - ^[5-^d{x, v)Y ) , 



bl{x) - i^l{x')]i^l{y)\ < ^ (^^) exp ( - j{S-'d{x, yDr ) exp ( - 7(<5-'d(y, vDT ) 

C /d{x,x')\v 



V{x,S'')\ (5'= 

and we may similarly sum over a € iV'^ . For S'^ < d{x, x') < d{x, y), we just use the quasi-triangle 
inequality and the first estimate of the lemma to both terms. □ 

9.2. Lemma. 

C 



< 



Vix,y)- 
Proof. By Lemma |9.1[ we have 



k,a k:S>'>d{x,y) ^ ' ' 

k:S>'<d(x.,y) ^ ' ' 

where the first part has the correct bound by Lemma 18.31 For the second part, we have 

k:5^ <d{x.y) 



< 



k:S'^<d{x,y) 

Coo ^ i 



E'^'"''exp(-7^— )<-— . □ 



9.3. Lemma. 



WAVELETS IN SPACES OF HOMOGENEOUS TYPE 



23 



Proof. By the second estimate of Lemma 19.11 

El[^'(^)-V'a(^')]V'S(2/)l 



k.a 



C fd{x,x')\'n 



k-.S'^ '>d{x.y) 



k:d{x ,x' )<.5^ <.d{x .y) 

The first two parts contain the factor x')/c?(x, y))'', and the rest is bounded by C/V{x,y) 
according to Lemma 18.31 and 'J2m=o (5"™*^*^"'"''^ exp(— 7(5^™°) < C. For the last term we even 
obtain the bound C/V{x, y) -exp ( — 7(o?(x, y)/d{x, x'))°- ) < C/V{x, y) ■ {d{x, x')/d{x, y))^ for any 
K. □ 

10. LiTTLEWOOD-PALEY decomposition and LP THEORY 

Recall that is the orthogonal projector onto and set also the orthogonal projector 
onto Vfe. These operators will provide us with a new regular Littlewood-Paley decomposition for 
spaces of homogeneous type. The following lemma describes the kernels of these operators: 

10.1. Lemma. The kernel Pk{x,y) of Pk is symmetric in x,y and has estimates 

y'^l{B{x,^'^))^I{B{y,^'^)) 



d{y,y')y ( cxp{~^i5~''dix,y)y) cxp ^i5~>^d{x,y')y 



\Pk{x,y) Pkix,y)\<C, J \^f,{B{xJ'^))f,{B{yJ^^)) ' ^f,iB{xJ'^))f,{B{y',S''))^ 

for some C, 7 and all x,y,y' G X and fc G Z (with k > ka if X is bounded). Moreover 

Pk{x,y)d^i{x) = 1. 

I X 

The kernel Qk{x^y) of Qk is symmetric in x,y and has similar estimates with s changed to a, the 
additional exponential factor 

exp(-7(J-'=d(x,^'^))'^), 

and the cancellation condition 

/ Qk{x,y)dfj.{x) = 0. 
Jx 

Note that there are many ways to express the denominator factors in the kernels up to chang- 
ing the constants C, 7 in the numerator factors especially thanks to the exponential decay. In 
particular, one possible expression shows that the system of operators Pk is an 'Approximation of 
The Identity' in the sense of |15[ Definition 2.2] (this part of that paper does not use the Reverse 
Doubling property), and the properties listed there hold. 

10.2. Theorem. The spline-wavelet representation yields a decomposition of Littlewood-Paley type 
with Holder- continuous kernels. 

Proof. If X is unbounded and / G L'^{X) the converging series 

00 00 
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is an homogeneous Littlewood-Paley decomposition. Also one can truncate at any level since 
Pk+i = Pk + Qk and write for any i 

oo 

k=l 

which is an inhomogeneous Littlewood-Paley decomposition. This decomposition is the one used 
if X is bounded with £ = ko. In that case, the term P^f is a constant. □ 

Observe that the sum J^k Qkf can also be rewritten as a "discrete" Littlewood-Paley decompo- 
sition J2k q(/i From there, one can look at convergence for / in various topological spaces, 
and develop the theory of function spaces. We restrict ourselves to spaces (and BMO in the 
next section) and leave further developments to the interested reader. 

The estimates in Section [S] immediately give the following result: 

10.3. Proposition. Let be arbitrary complex coefficients bounded in absolute value by one. 
Then the series 

fceZaeAfc 

converges absolutely for x ^ y and satisfies 

\Kix,y)\ < \Kix,y) - K{x',y)\ < ^-(^p^Y , d{x,x') < 

V[x, y) V{x, y) V d{x, y) J 2Aq 

with a similar Holder- regularity estimate in the second variable. 

10.4. Corollary. The spline wavelets form an unconditional basis of L^{ii) spaces when 1 < p < oo. 

Proof. Completeness follows using the convergence properties of the Pk in LP{ii). It remains to 
show that operators Tc given Tc{tp'^) = c^ip'^ are uniformly bounded on all LP{h) spaces whenever 
c = (c^) is a sequence of complex numbers in the unit ball of C. These operators are contractions 
in L^ifJ.). Using the regularity of their kernels as proved in Proposition 110.31 and the Calderon- 
Zygmund theorem, Tc has weak type (1,1), with uniform bound with respect to c. Since T* = Tg, 
the same applies to T* . We conclude by interpolation. □ 

11. BMO THEORY 

Recall that the space BMO(/^) of functions of bounded mean oscillation consists of those 
b^LlM with 



I^IIbmo(m) ■= supinf — ^ f \b-c\d^i< oo, 



where the infimum is almost realized by c = &b := iJ,{B)^^ jgbdfi. We do not incorporate 
the norm of constants as a part of the BMO norm even if X is bounded, so that our BMO is a 
Banach space modulo constants in both bounded and unbounded cases. These averages satisfy 
the following useful estimate: 

11.1. Lemma. Let Bi = B[xi,ri), i = 1,2, be two balls in X. Then 

ri +r2 + d{xi,X2) 



\bBi - 1 < II&IIbmo(m) (l + log ■ 



min{ri,r2} 



Proof. Suppose first that ri < r2 and Bi C B2. Let := Bi and be the smallest 2^ B"^ such 
that /i(S*+i) > Cf,fi{B^); hence also ix{B^+^) < Cf,fi{B'). It is easy to check that - < 
||6||bmo(ai)- Let B'' be the first B^ with radius bigger than that of B2. Then we also have 
\bBk - ^ II&I1bmo(m)' aiid A: < 1 + logr2/ri. 

In the general case, we may choose an auxiliary ball B^ of radius r^ ^ ri +r2 + d{xi, X2) which 
contains both Bi and B2, and apply the earlier consideration to 
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11.2. Corollary. Let f be a function with exponential decay \f{x)\ < Ccxp(— c • (i(a;, Xq)") with 
C,c,a > 0, and b G BMO(/i). Then the product f ■ b is integrable over X. 

Proof. Let Bn B{xq, n). Clearly /. and hence f -bs^ is integrable, so we consider f ■ {b~ bsi). 
By Lemma ril.li 

/ \b-bBi\d^i< \b-bBjd^ + fi{Bn)\bB„-bBi\<f^{Bn)il + logn)\\b\\BMO(p), 

J Br, JB„ 

and li{Bn) ^ n*^/i(i?o) by doubling. Hence \\lB„{b— ^So)lliHAi) grows at most polynomially in n, 
whereas ||1s„\s„_i/|1l=c(^) decays exponentially in n. Thus f ■ {b — ^So) ^ ^^(A')- ^ 

This implies in particular that the wavelet coefficients {b.-tp'^) = J-^ b ■ ip^dfi are well-defined 
for b G BMO(/i). The following injectivity property is somewhat technical, and we postpone its 
proof to Appendix [X] 

11.3. Proposition. Suppose that b E BMO(/^) satisfies (&, V'q) = for all k E Z and a E '3/^^ . 
Then b is equal to a constant. 

If X is bounded, one should mention that k > ko, but in fact, notice that ^^'^ = and so ijj^ 
does not even exist when k < ko. We thus do not need to distinguish further between X bounded 
or not. 

Note that some nontrivial a priori size condition on b is in general necessary for such a conclusion. 
For example, if the ip^ are regular wavelets on K*^, then ('0^, P) = for all polynomials P of degree 
lower than the regularity of the wavelets. 

We say that a sequence {^alfcGZ.ae®''' ^ Carleson sequence if 

||{6^}fc,a||car := SUp (— ^ J2 I I ') < ■ 

(fc+l,a)<(«,/3) 

Pay attention to the fact that the supremum runs over all £ and f3 E which index the dyadic 
cubes Q^, whereas the sum runs over k and a E = . Via the wavelet decomposition, 

we obtain an isomorphism between BMO functions and Carleson sequences: 

11.4. Theorem. The spaces BMO(/^)/C (13M0 functions modulo constants) and Car are isomor- 
phic. This isomorphism is realized via b h- !■ {{b,ip^)}k,a, with inverse given by 

(11-5) {bt}k,a ^ ^"(^' - hr.s^>ro}ik)^'M) , 

k.a 

where the series converges in L^^^^n) for every xq E X and tq > 0, and the choices of xq and rg 
ordy alter the result by an additive constant. 

The result and its proof are reasonably classical in spirit, but the lack of bounded support of 
the wavelets somewhat complicates the matters. While exponential decay is intuitively almost as 
good, one needs to go through certain technicalities if one wants to be careful with convergence 
issues. We indicate the argument. 

Note that we have formulated the recovery of the BMO function from its wavelet coefficients 
using the "infra-red" renormalisation of the wavelet series (since the modification on the series 
appears for large scales which corresponds to small frequencies in the classical Euclidean situation) 
rather than the iJ^-BMO duality; this is in contrast to the statement in Meyer's book [30], for 
instance, but this other method would be certainly doable here as well. 

Proof. Let us first see that b ^ {{b^^a)}k,a maps BMO(/x)/C into Carleson sequences. The 
injectivity of this mapping is the content of Proposition [TT731 Given (fc, a), we let B ~ B{x'^, CS^), 
with some large C, be a ball such that d{Q'^, (B^Y) > 6''. We write 

b={b- bg)l B + {b- bs)lB^ + b^+b^ + b\ 

Then 

Y.\{b\^pf,)\^<\\bX<\MBMOMKQi) 
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by orthogonality. John-Nirenberg inequality and doubling; (&^, ip^) may be estimated by the decay 
of the wavelets, and {b^,ip'^) = 0. 

From Carleson sequences to BMO. Conversely, assume the Carleson condition. Given a ball 
Bi = B{xi,ri), we may rearrange 1^^ times the right side of (jll.Sp as 

E E &'^'' + ib, e E bi^a + ^B, E E^'(^'-^'(^i)) 
-ib, e Y.bai^'M + iB, E E^'^'(^i) 

k:ro<S''<ri a k:ri<5''<ro " 

+ Ib, E &a(V'^(2;i)-V'S(a:o)). 

^.i5fc >ro Vri 

Here the last three terms converge uniformly to Ib^ times a constant, the first term converges 
in L'^in) with norm bounded by (use the Carleson condition after covering CBi by 

boundedly many dyadic cubes of sidelength ~ ri), and the second and third terms can be estimated 
uniformly by the decay and regularity of the wavelets ip^. This proves the L^{fJ.) convergence on 
Bi, and also the BMO estimate 

||lB,(&-CBj||2<M(i?l)'/', 

where b stands for the function on the right of pi.Sp . and cbi is the constant produced by the 
last three terms in the above expansion. Effectively, the same argument also shows the possibility 
of replacing (xo,ro) in (|11.5p by (xi,ri), only changing the result by a constant. 

Verifying that the two mappings are inverses to each other. This is the most technical part of the 
argument. For a Carleson sequence {b^}k.a, let b denote the BMO(yL/,) function (as shown in the 
previous part) on the right of (jll.Sp . We claim that {b,ijj^) = 6^, which completes the proof. 
Indeed, this also shows that if {b^}k,a arose as 6^ = (bjip^) from some b E BMO(/x), then the 
function b produced by (jll.Sp satisfied {b — b,ip^) =0 and hence b — b + constant. 
While the claim is formally obvious, due to the orthogonality 

we need to justify exchanging the order of summation and integration. To this end, let us reinves- 
tigate the convergence of (|11.5p where, since -0^ is orthogonal to constants, we may assume that 
(xo,ro) = iy^^,S''). Then 

^ = E E E Y.b'M-^Uy'))=--b' + b'- 

{i+i,i3)<{k+i,e) e:S<^>s'' p 

The second part converges pointwise absolutely to a limit of size < 1 + \ogj^{d{x , y%) / S'') . Since 
■0^ decays exponentially, this allows to exchange the summation and integration when computing 
(6^,-0^). Since is orthogonal to for 5^ > S'^, as well as to constants, we get (&^, V'a) = 0. 
For the first part, we have 

61 =: E = E E Uq-^^MU < y^^Op) exp ( - 7('5-'d(4-+\4+'))'^ ), 

e e c '' '' 

whereas 

l|l0j+^V'al|2 < CXp ( - 7(rt(4-+\x^l)r ). 
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Hence 

e e c 

< E E \/M(QF^exp ( - 7(<5-'=d(4+\ 4^+1))'^ ) exp ( - x^^))" ) 

;S E exp ( - ^{6-'^d{x'^+\xt+')r ) < 

C 

This allows the first exchange in 

(&\V'^) = E(^'^'^'') = E E ^^('/'^-V'^) = E E bise,s^p = bi, 

where the second exchange follows from the convergence of fi)<(k+i 9) ^^V"^ in -^^(m)i a-^d 

from -0^ S ^^(m)- Altogether, we get {b,il!^) = 6^, as claimed. □ 

12. The r(l) theorem 

To illustrate the power of the spline wavelets, we use them to sketch a proof of the T(l) theorem 
in any space of homogeneous type. Such a result is surely part of the folklore, but surprisingly 
difficult to find spelled out in complete generality: the seminal paper of David-Journe-Semmes 
[S] makes several assumptions on the space, like the small boundary property of balls, and many 
recent references treat other special cases like Ahlfors-David [TT| or reverse doubling spaces [15] . 

The technology to prove the T{1) theorem in a general space of homogeneous type has certainly 
existed since the work of M. Christ Indeed, a proof of the T(l) theorem can be given by using 
the Haar wavelets only, and these have been available since Christ's construction of his dyadic 
cubes with the small boundary property. In fact, Christ even formulates the general r(l) theorem 
[S] Theorem 8], but attributes it to [H], and proceeds to use it as a tool for proving a certain 'local' 
variant. However, Christ's techniques would have clearly delivered a proof of the 'global' T(l) 
theorem as well, without the restrictions imposed in [9]. 

So the regular wavelets are not strictly necessary for obtaining the T(l) theorem, but they 
nevertheless provide a rather efficient tool for that purpose. We only indicate the argument, 
which largely imitates the treatment in Euclidean cases given by Meyer [3D] . 

We take as space of test functions the space = Cq{X) of functions with bounded support 
and Holder- regularity s equipped with the usual topology, where s € (O,?^) is arbitrary and we 
recall that 1] is the regularity of the splines. The space Vs is dense in by Proposition 14.51 

Let V'g denotes its dual space. Recall the standard definition. 

12.1. Definition. Let T : be a linear continuous operator. We say that T is associated 

to a Calderon-Zygmund kernel of order s if the distributional kernel K{x, y) of T satisfies for some 
constant Ci < oo, 

\K{x,y)\<CiVix,yr^ 

when X ^ y and 

\K{x,y) - K{x,y')\ < Ci (^^^^J V{x,y)-^ 
when x,y,y' S X with < d{x,x') < dix,y)/{2Ao) 

\K{x, y) - K{x, y')\ < C, f ' V{x, yY^ 

when x,x',y £ X with < d(x,x') < d{x,y) / {2Aq) and if furthermore, for any / e Vs, one has 
the representation 

(12.2) Tf{x)^ I K{x,y)f{y)dpi{y) 

for almost every x ^ supp/. 
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Recall also the weak boundedness property and the meaning of r(l). A linear continuous 
operator T : Vg ^ has weak boundedness property WBP(cr) if 

\{Tf,g)\<CoV{x,r) 

whenever f,g € Va, with support in B{x,r) and are normalized by ||/||cxj + '"'^ll/llcf ^ 1 a-nd 
similarly for g. It is well-known that WBP{(t) and WBP{<t') are equivalent conditions whenever 
s < (t,<t' < 7] when T is associated to a Calderon-Zygmund kernel of order s. As for T(l), it is 
defined as a continuous linear functional on the subspace of of functions / with mean value 

by 

(T(l), /) = {Tg, f) + jj^- 9{^)) CTf){x) dfi{x) 

with g a function in Vs that is 1 on a ball B(xo,r) containing the support of / and on the 
complement of B{xo, 2Aor) as in Corollarv l4.2[ as *T/ is integrable away the support of /. 

12.3. Theorem. Let (X, d, /i) be any space of homogeneous type andT be associated to a Calderon- 
Zygmund kernel of order s. Then T has a bounded extension to L'^ilJ.) if and only ifT has WBP{s), 

r(i) G BMO(Ai),*r(i) e bmo{^i). 

As usual it suffices to prove the converse. With our spline wavelets at hand, any of the stan- 
dard wavelet proofs of the T(l) theorem applies to our statement. For example, one can follow 
almost line by line the wavelet proof given in \50\ pp. 267-270]. The weak boundedness prop- 
erty and the kernel representation allows to define {T{l),ip^), whose absolute value has a bound 
Cy/ ^{B{y^, 6'^)). Then, one takes away the paraproducts and set U = T — IiT{i) — ^^tT{i) where 

with s^+^ the spline s^"*"^ normalized in Recall that is localized near the point y'^ = 

j.k+1 g a^k _ (^k+i \^ i^k^ gQ ^j^^^ ^]^g Calderon-Zygmund estimates follow from the results in 

Section |9] (see the initial comment there) when |(/3,'0q)I 1^ C-\/ fi{B{y^, (5'"')). The paraproduct 
is classically bounded on L^{n) if and only if the coefficients form a Carleson sequence. 

It uses the maximal function f*{x) = sup |(/, s^+^)| taken over all k G Z,!^"*"^ € such 
that X £ B{x'^^^ , CS'^^^) which is clearly comparable to Hardy-Littlewood maximal function. By 
Theorem lll.41 this is equivalent to /3 S BMO(/x)/C. Thus U is associated to a Calderon-Zygmund 
kernel of order s with, moreover, U{1) = *f/(l) = 0. 

Next, by taking /, g as finite linear combination of the wavelets (which are dense), it suffices 
to prove that the coefficients of U on the wavelet basis form a bounded matrix on Here, using 
only kernel regularity estimates (not the size, in fact), weak boundedness, U{1) = *C/(1) = 0, one 
makes sense of the coefficients (Uip^jipp) and finds the estimate 

(12.4) |(f/V„>^)| < |,iB{y^^,S'^)) + ^,iB{yf,,^^)) + V{y^^,yf,) 

with < e < s and k A £ ~ inf(fc,^). The version of the Schur lemma for the i'^ boundedness is 

E |(C/^^,4)lv/MS(yS,'5'^-)) < ^ti{Biyl,S^)) 

{k,a) 

uniformly in {£,/3) and the symmetric estimate reversing the roles of {k,a), {£,(3). Details are left 
to the reader. 

Strictly speaking, this argument works when n{X) — oo. When fi{X) < oo, one has to incor- 
porate the following observation. First, the assumption that T(l) G BMO(^) implies in particular 
that T(l) is (identified to) a locally integrable function. Second, since X is bounded, it can be 
regarded as a ball, and the constant function 1 is a smooth bump function associated with this 
ball. Hence, the weak boundedness property implies that | J-^ T{1) d^\ ~ |(r(l), 1)| < \i.{X\ and 
thus |1T(1)||bmo(ai) + I /x ^(-'^) "^/^l < °° (recall that our BMO norm is the homogeneous norm). 
By John-Nirenberg's inequality, this implies that r(l) e L'^{\s) with ||r(l)||2 < ||T(1)||bmo(m) + 
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\J^T{l)dfi\. Similarly ||*r(l)||2 < ||'r(l)||BMO(M) + 1 /x ^A^l < °°- The argument before can 
be repeated and implies that ttTtt is bounded where tt is the orthogonal projection onto the sub- 
space of functions in with mean value since the wavelets span this space. The boundedness 
of T on (p) follows readily. 

12.5. Remark. As in |30| . estimate (|12.4p is stable under matrix multiplication up to changing e to 
a smaller value. This algebra property furnishes a proof that Calderon-Zygmund operators with 
T(l) = *T(1) = is an algebra for the composition. This was proved in [T3j by working with a 
new Holder-continuous quasi-distance as in [5], hence changing the class of singular integrals as 
discussed in the Introduction. 

13. Redundancy of the size estimate 

It turns out that the size condition on the kernel is actually redundant in the T(l) theorem, 
in that it already follows from regularity of the kernel and the weak boundedness property. To 
our knowledge, this remark seems new even in the context of 'R'^ with the Lebesgue measure. We 
came across this observation noticing that the regularity estimate of kernels like J2 ^a{^)^a(.y) 
with /i^ satisfying the size and regularity of spline or spline- wavelets is straightforward, while the 
size estimate required the analysis in Section IH] based on the structure of the point sets . It is 
also possible to develop the T(l) theory without using the size estimate at all. 

To be more precise, let us say that T : Vg ^ Vg is associated to a Calderon-Zygmund kernel 
of order s in the relaxed sense, if in Definition 112.11 we leave out the condition that |/4r(a;,y)| < 
CiV{x,y)~^ , and only assume that K{x,y) is locally integrable away from the diagonal. This 
condition could be further relaxed, by not assuming the a priori existence of a measurable kernel 
at all, only that T is a weak limit of operators T„ . which are associated to relaxed kernels of order 
s in a uniform way. However, we stick to the stated relaxation, for the simplicity of formulation. 

13.1. Proposition. Suppose that T : Vs — > is associated to a Calderon-Zygmund kernel of 
order s in the relaxed sense, and that T satisfies WBP{a) for some a G [s, 77]. Then T is actually 
associated to a Calderon-Zygmund kernel of order s in the sense of Definition \12TT[ 

Proof. Recall that the paraproducts n7-(i) and *Ht7-(i) are associated to Calderon-Zygmund kernels 
of order s under these assumptions. Thus it suffices to prove the claim for the operator U ~ 
T — \It{i) ~ *nt7-(i) in place of T. The kernel of U is given by 

K{x,y)^ ^ {U^l4)^t{y)4i^l 

k,i,a,P 

where the coefficients satisfy the estimate (|12.4[) , as this estimate did not make use of the size of 
the kernel. 

We show that this series converges absolutely and satisfies the size estimate. By symmetry, it 
suffices to consider the half of the sum with £ < k, hence > 5'^, and 

\{Ui;l4)^l;'^{y)4{x)\ 



< 



S(>^-^>il + 6-'d{ytx$))- ^ ^ 



■ exp ( - j{5-'^d{y, ) exp ( - i{5-H{x, y^T ) 



Hence 

5^ Y,m^l4)^'^iy)4{x)\ 



k,t.t.<k a„ 



i:S'>d{x,y) k:S''<S' ^ ' k,e:d{x,y}>S'>6k ^ 
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where {l + 6^^d(x, ?/))"'* was dominated by 1 in the first sum, and by {S~^d{x, y))^" — S^'^d{x, y)~^ 
in the second. In the first part, we sum a geometric series in k, and then use the bound of 
Lemma 18.31 In the second part, we simply sum up a geometric series in k producing the factor 
iJ^", and then a geometric series in ^. □ 



Appendix A. Proof of Proposition 111.31 
We begin with a technical estimate: 

A.l. Lemma. For a any ball B — B{xo,r) and b G BMO(/.i) of norm one, we have 

1 + log (((5'' + r + d{y^,xo))/ mm{S'',r}) always, 

2Tq 



\^ iJa^ 2A, 



B and S'' < r. 



with B^ :=B(y^,(5fc). 



Proof. By estimating l-iA^I < C/i(B^)"i/2g-7'5 \ S-^"'-^^B^, and simply ignoring 

the indicator, we estimate up by 

{\^':^\,\b-bB\)<^^Y.e-^'-"'^ [ |6-5B|d/i, 
\/u(B^) Js-^^Bi 



where 



/ \b-bB\dfl< f \b~bs-r.B>'\dfl + fi{S-"'B^^)\bs-r.B^-bB\ 

V mmjd ™o'',r)| 

5'' + r + d{y^, xq) 



<r"^/i(S^)(l + m + log 

and hence 



(|V'SU6-6B|)<C^M(i3S)(l + log 



min{(5'^, r} 

S'' + r + d{y^,xo) 
min{6^, r} 



If y^ e JJ^-^' ^^^^^ d{y^, B^) ^ r, and 5^ < r, we make a similar estimate but exploiting the 
fact that for S-"'B^ to touch B", we need S-'"'d'' > r{B). Thus 

i\i^'^llB^-\b~bB\)<^^ e--''-"'^ f \b-bB\d^,, 



^U§m ^ e---M.-™i?^)(l + ™ + log^ 



7n:S >r / 



We are ready for the proof of Proposition 111.31 which we recall here: 

A. 2. Proposition. Suppose that b E BMO(/i) satisfies b) ~ for all k, a. Then b is equal to 
a constant. 

Proof. We show for every ball B = B{xo, r) and e > the following: there exists a constant c such 
that ||ls(6 — c)||oo < This clearly suffices. We may assume for simplicity that ||6||bmo(^) < 1- 
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Given B, we take a large auxiliary B = i?(.To,f). We use the fact that lg{b — bg) E can 
be expanded in terms of the wavelets ip^: 

k.a 

(A.3) =is E E^'W'- 15(^-^5)) + Is E E(^'-'^'(^o))W'-is(^-^s)) 

k:S^<r' a k:S''>r' 01 

+ is E E^'(^o)('^-ii3(^-^B)) + is&B. 

Provided that everything converges (which we check in a moment), the terms on the last line give 
1b times a constant, so it suffices to show that the second-to-last line becomes arbitrarily small 
for properly chosen f > r' > r. For the convergence of the last line, note that 

\^tixo)i^llB{b^b^))\ < mXo)\Ut\U\lB{b-b^)\\, 

c 



< 



and, by Lemma 



exp(-7(<5-*^d(xo,y^))'')MS), 



E E77§iT^-P(-7(rt(xo,2/S)r) 

So the last line of (jA.3[) is a well-defined constant, as claimed, and it remains to see that the rest 
of the right side of the same equation is small. 



Part 6'' < r' . Since V'q is orthogonal to constants always, and to b by assumption, we have 

{i^iAB{b-b^))^~{i,li^.ib-bs)), 

and we may apply Lemma [A. 1 1 f with B in place of B) to estimate. If G {^^^BYi then 



UllBAb^b^))\<^jli{Bi){l + \oi 



d{ya,xo] 

Sk 



whereas 

II Ib^^IIoo < exp ( - j{6-''dixo,y^)r ) ■ 

V Ml-Daj 

Hence 
(A.4) 

5^ ^ \\lB^t{^PtlB{b-b^))\U< cxp(-7(rt(^o,2/^)r) 

k:S''<r' 

If, on the other hand, g 2J^B, then we have the estimates 



\ii^llB^ib-bB))\<^/M)-e-^^'/''^\ 



< 



hence 

E E \\^B^'MiAB{b-b^))u< y: E 

(A 5) ^'■^''^''''^■■v!i^2h^ ''■■■^''<^'<^--ylieji^B 
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Part S'^ > r' . As before, we have from Lemma [A. II (with B in place of B) that 

isib - b,))\ < v^(i + log +; 

+ iug+ — -I- ^"8+ 

r 



< Vm(SS) ( 1 + log+ - + log+ ^ + log+ 



and from the regularity of the wavelets that 

II - i^'M)\\oo < (^) ' exp ( ~ 7(<5-'rf(?/S, ^^o))"^ ) . 



Hence 



k:S''>r' a 

r \v / f 5^ 



(A.6) 



fe:i5'' >r' 



X ^(1 + log^ *|lf2l)cxp( - 7(i-'<l(»J,io))"). 

Conclusion. Substituing (jA.4p . (jA.SP and (|A.6P into (jA.3|) . we have shown that: for any B — 
B{xo, r), and r > r' > r, there exists a constant c = c(_B, f , r') such that 

||lB(&-c)|U<e--('VO" + (^)''(l + log^). 

If we now choose r' = tr, f ^ tr' = t^r, we get 

||lB(6-cs,t)l|oo<e-''*"+r''(l + logt) 0. 

t— f oo 

This concludes the proof. □ 
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